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A branching random walk in presence ol an absorbing wall moving at a constant velocity v 
undergoes a phase transition as the velocity v of the wall varies. Below the critical velocity Vc, 
the population has a non-zero survival probability and when the population survives its size grows 
exponentially. We investigate the histories of the population conditioned on having a single survivor 
at some final time T. We study the quasi-stationary regime for v < Vc when T is large. To do so, one 
can construct a modified stochastic process which is equivalent to the original process conditioned 
on having a single survivor at final time T. We then use this construction to show that the properties 
of the quasi-stationary regime are universal when « — > Uc. We also solve exactly a simple version 
of the problem, the exponential model, for which the study of the quasi-stationary regime can be 
reduced to the analysis of a single one-dimensional map. 



I. INTRODUCTION 



Branching random walks are often used as simple models of evolving populations, with or without selection. They 
can describe how a population invades an empty domain, how a favorable mutation spreads or how the fitness 

of a population evolves under selection [l^, [23 . Recently very simple models 0, H, [l^, El of neo-darwinian evolution 
[3^ have been studied, motivated in particular by in vitro experiments p^ . In these models each individual of a 
population of fixed size is characterized by a single number, a trait which we can call its fitness or its adaptability. 
From one generation to the next, this trait is transmitted to the offspring up to small variations which represent the 
effect of mutations. Then at each generation, only the N fittest individuals survive. It was shown in that these 
models can be mathematically reduced to the study of travelling waves in presence of noise, a subject widely studied 
[T^ . [TtI . [sp , m, [s^, [37I. 1 3 SI which appears in many other contexts such as reaction-diffusion problems, disordered 
systems Q or QCD [H, M [11 . 

Another version of the same model of evolution was also considered in [l^ where the size of the population may 
vary with time and the selection takes place by eliminating at each generation all the individuals which are below 
some threshold along the fitness axis. When this threshold increases linearly with time, the model reduces to the 
problem of a branching random walk in presence of an absorbing wall moving at a constant velocity v. Depending on 
how fast the wall moves, the population can either survive or get extinct. The survival probability, in the long time 
limit, decreases as the velocity v of the wall increases and, above some critical velocity Vc, the population gets extinct 
with probability 1. 

A model fully equivalent to the branching random walk in presence of an absorbing moving wall, that we study 
here, is the case of a branching random walk with a drift and a fixed absorbing threshold, which has been introduced 
in [l[ in order to describe the aging process in a cell proliferation model. In this case too, there is a phase transition 
to an absorbing state when the drift increases. 

Similar phase transitions to an absorbing state (the empty population) have been extensively studied in the th eory 
of non-equilibrium systems (2^ . [36j . in problems like directed percolation [15] or reaction-diffusion processes [26| . 
Beyond the characterization of the transition to this absorbing state, there has been an increasing interest to study 
the process conditioned on the survival at a very late time [ll|, [ij, [24I . 

In many cases, evolutions conditioned to avoid the absorbing empty state lead to a quasi-stationary regime : for 
example, in a birth-death process (see [H m, 113, HI and section |TT] below), if one conditions on a given non-zero 
size at a final time T and looks at the population at a time t such that both t and T — t are much larger than the 
characteristic relaxation time of the system, the statistical properties depend neither on the initial condition nor on t or 
T and one observes a quasi-stationary state^ . The probability of the events which contribute to these quasi-stationary 
regimes decreases with T, and therefore they are difficult to observe in simulations. There are however a number of 
cases for which one can modify the process (i.e. determine the right biais) [2l|, to generate the quasi-stationary 
regime directly. 



In the present paper, we use the term "quasi-stationary" to describe a situation at time t conditioned on the fact that there is a fixed 
non-zero number of individuals at a much larger time T. Note that "quasi-stationary" is often used to describe a different situation 
where the two times t and T coincide |13i . 1I81 . i45f] . What we call here a quasi-stationary process is sometimes called a "Q-process" [lO| . 



Figure 1: Branching random walk in presence of an absorbing wall. The circles represent branching events and the crosses 
absorption by the wall. 



In the present work, we consider the problem of a branching random walk in presence of an absorbing boundary 
which moves at a constant velocity v (see fig. [T] and [13] )■ On the infinite line, i.e. without any boundary, the 
population spreads linearly in time [Slj at a velocity Vc and the number of walkers grows exponentially with time. 
The presence of a moving absorbing wall introduces a competition between the growth and the extinction due to the 
wall. When v > Vc, the population gets extinct with probability 1, whereas when v < Vc, this extinction probability 
is strictly less than 1 and there is a non-zero probability that the population grows exponentially with time. Here, 
we mostly study the quasi-stationary regime, conditioned on the survival of a few individuals at a very late time T, 
for V < Vc- This quasi-stationary regime has already been studied in the mathematical literature, where a modified 
process was invented (2ll [22] which allows one to generate typical evolutions in this quasi-stationary regime. Our 
main goal in the present paper is to show how the properties (average size of the population, density profiles) can be 
calculated in the quasi-stationary regime and how universal expressions for these properties emerge when v ^ Vc- 

The paper is organized as follows : in section |TT1 we recall how the quasi-stationary regime can be calculated in 
the case of a simple birth-death process, and how the distribution of the population size becomes universal near the 
transition. In section [Till we extend the approach of section [TTj to the branching random walk process in presence of 
an absorbing wall. In section IVl we show how the quasi-stationary regime can be generated by a modified process, 
as already described in [U, [22| and we obtain explicit expressions of the density in this regime. In section |Vl we show 
that the average size of the population and the average density profile become universal as v — + v^. In section [VTl 
we analyze a simpler model, the exponential model, for which the whole distribution of the population size can be 
calculated. 



II. AN EXAMPLE OF BIRTH-DEATH PROCESSES : THE GALTON- WATSON PROCESS 

In this section we discuss the simple case of the Galton- Watson process [1^, [l^l . Our goal is to explain in this 
well-known example the approach that we will use in a more general context in sections IIIIl IIV[ FVl 



A. Discrete time 



We first consider the discrete time case. The population is fully characterized by its size Nt at time t. At every 
time step, each individual is replaced by k offspring with probability pk (in particular, po is the probability that an 
individual leaves no descendance). Let us define the extinction probability Qe{i) as the probability that iV; = 0, given 
that we start with a single initial individual at time 0. During the first time step, the initial individual may branch 
into k offspring : thus, extinction of the initial lineage after t + 1 time steps is equivalent to the extinction of the 
lineages of its k offspring after t time steps. Since in this model the descendants of these offspring are independent, 
Qe{t) satisfies the following exact recursion : 



Qe{t+l)=F{Q,{t)) 



(II.l) 
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Figure 2: Size Nt of the population at an intermediate time t knowing that the final size is Nt = m at time T = t + t' . 



where 

oo 

-ifc 



FiQ)^yPkQ' {11.2) 



with the initial condition Qe{0) = 0. For large t, Qe{t) converges to a limit Ql and 1 — Ql is the probability of eternal 
survival of the population. This limit Q* is the attractive fixed point of the map F, satisfying Q* = F{Q*). The 
function F{Q) is convex and satisfy F{\) = 1 and F{0) = po > 0. The average number of offspring of an individual 
is given by F'{1) = J^kPk^ — k. If fc < 1 the stable fixed point is Q* = 1 and the population dies with probability 1. 
On the other hand when fc > 1, there is a non-zero survival probability 1 — Q* in the t ^ oo limit. 

If one starts at i = with a single individual, the distribution of the population size Nt can be characterized by its 
generating function 

Gi(t;Ai) = (e-''^'). (II.3) 

The size Nt is the sum of the contributions of all the offspring dX t — 1. The independence of the lineages of the fc 
offspring implies that Gi evolves according to the same equation (jll.l[) as Qe{t) : 

Gi{t+l-ti)=F{Gi{t-ii)), (II.4) 

the only difference being in the initial condition Gi(0; /i) — e^^. Note that the extinction probability is given by the 
particular case Qe{t) — Gi(t;oo). 

For large T and a given finite size m, the events Nt = m are more and more rare since the population size 
either vanishes or grows exponentially. For such events, however, one may be interested in the sizes Nt, < t < 
conditioned on the fact that Nt = m (as in figured]). To do so, let us divide the time interval [0, T] into two intervals 
of lengths t and t' = T — t and consider the two-time generating functions G2{t, t'\ fi, v) defined as : 

G2(t,t';M,i') = (e-'''^'-'''^'+">- (11.5) 

As for Qe and Gi, one can easily show, by considering the first time step, that G2{t,t'; fi, v) as a function of t also 
evolves according to pi.ip : 

G2{t + 1, t'- fi, ly) = F{G2{t, t'; fi, i/)), (II.6) 

with the initial condition 

G2(0,t';Ai,i.) = e-^Gi(i';i^). (II.7) 
If one expands G2{t^ t'] ^, u) in powers of : 



oo 

E 

m=0 



G2{t,t'-iJi,v) = V e-''™i?„,(i,i';Ai) (IF 



the coefficients Rm{t,t'; fi) are related to the generating functions of Nt conditioned on the survival of m individuals 
at time T (i.e. Nt = m) : 

Rmit,t'; ^i) = P^it + t'){e-''^'\size Nt+v = m) (n.9) 
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where Pm{t + t') is the probability that iVf+f = m. As in the special case /i = 

oo 

G2{t,t';0,i^)^Gi{t + t';i^):^ ^ e-'^"P„,(i + i') (11.10) 

one can extract the probability Pm{t + t') from the expansion (jll.lOp of Gi{t + t', v) in order to normalize Rm{t,t'; /i) 
in (jll.Qp and get the conditional probability (e~^^*|size m a,t t + 1' ). 

We are now going to see that (e"''^* |size m at t + t' ) has a limit as t and t' ^ oo and to show how to compute 
the distribution of in this quasi-stationary regime. Since Gi and G2 evolve according to (jII.4III.6|) , their long time 
behaviours can be extracted from the properties of F near the fixed point Q* . If ut is a sequence obtained by iterating 
a map F 

ut+i = F{ut) (11.11) 
with an initial condition uq. If it converges exponentially to a fixed point Q*, one can write : 

wt ~ Q: + AU(uo) +o(A*), with A = F'(g:) (11.12) 
where the amplitude A{uq) is defined by : 

^(«o)- lim (11.13) 

Using (|II.12[) for G2 and Gi together with the initial condition (|II.7p for large t and t' gives : 

G2{t,t'-n,v) = Q:+A*y^(e-^Gi(i») + o(A*) 

= 0: + A*yl (e-^Q: + e-^K'' Aie-") + ...)+... (11.14) 

= Q: + K'A (e-^Q:) + A*+*U' (e-'^Q:) e-^yt(e-'') + . . . 

If one expands ^(e^") = A.^e"'^"^, one gets from (111.81 III.10III.14|) at leading order in A*+* for m > 1 : 

P^{t + t') ~ A*+*U'(0:)A„ = A*+*'A™, (IL15a) 
Rra{t,t'\lJ.) ^ A*+*'^'(e-'^Q:)e-''A„. (II.15b) 

From these equations and pi.Op . it is easy to see that for t,t' ^ 00 : 

(e-'^^'lsize Nt+t, = m ) (e-^'^)q. 

with 

(6-"^)^, = ^'(e-''g:)e-^ (11.16) 

This shows that the knowledge of the amplitude A defined in (|II.13p determines the whole distribution of the population 
size in the quasi-stationary regime (28j . 

In general there is not a closed expression of A for a given map F. One can however relate the expansions of A 
and of F in the neighbourhood of the fixed point Q*. The definition pi.lSp implies that A satisfies : 

AAiQ) = AiF{Q)), (11.17) 

a{q:) - 0, a'{q:) = 1. 

By taking successive derivatives of pi.l7p at x = Q* , one can relate the Taylor expansions of A and F 

-4(q: +6) = 6+^^(2) + + 

2 3 

F{Q: + e)^Q:+Ae+ '-F^^^ + '-F^^^ + ... 
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and get : 

^(^) = (IL18a) 

(3) P'-^^ 3(F(2))2 

" A(l-A)(l + A) + A(l + A)(l-A)2 ^"-^^^^ 

etc. Then, the moments of Nt in the quasi-stationary regime can be obtained by taking successive derivatives of 
(|IlT6|) at = : 

(iV)q, = 1 + Q:^(2) (n.l9a) 

{N^)qs - 1 + 3Q:^(2^ + (Q:)M(3\etc. (II.19b) 

As one approaches the transition point (i.e. k ^ 1 and Q* — > 1), expressions ()II.18|II.19P diverge since A ^ 1 and 
the moments scale as : 



F{2) 

(7V2)q, ^(3) 3A (1~A)F(3) 



(11.21) 



{N)ls (^(2))2 (i + A) (i + A)(F(2))' 

and so on for higher moments {NP)qs/ {N)p^^, p > 2. When the first moments of the number of offspring are finite (for 
example when k'^,k^ < oo), the ratios such as (1 — K)F^'^^ / {{1 + A)(F(2^)^) go to zero as A ^ 1 and one gets : 



(iV) ~ ^^£Z^ (II.22a) 

^ l-F'{Q*) ^ ^ 

^^'^'i^ • ^ (II.22b) 



{N)ls 2 

3, etc. (IL22c) 



(N) 
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One can notice that these expressions are valid on both sides of the transition, i.e. for fc — ^ 1+ and k ^ 1~ but do 
not depend otherwise on the Pfc's. This agrees with the fact (see [2^ and section Hi CI below*) that the distribution of 
Nt in the quasi-stationary regime becomes universal near the transition k — 1. 



B. Continuous time version 



We now consider the continuous time version of the previous birth-death process : at any instant t, the population 
is fully characterized by its size Nt. During every infinitesimal time dt, each individual either dies with probability 
f3odt or branches into k offspring with probability Pkdt. If one divides the time interval [0,t + dt\ into the first interval 
dt and the remaining interval [dt,t + dt], the independence of the evolution of the offspring implies that the extinction 
probability Qe{t) satisfies as in (|II.1|) : 

Qe{t + dt)= 1 - 51 l^kdt Qe{t) + ^ l3kdtQe{tt (11.23) 

y fe>o I k>o 

which gives in the limit — > : 

dtQeit) = F{Qe{t)), F{Q) - /^i^iQ^ - Q)- (11-24) 

fc>0 

This is the continuous time version of (jll.ip . In the long time limit, Qe{t) converges to the stable fixed point Q* of 
F : 



FiQ:) = 0, F'(g:) = A<0. 



(11.25) 
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The transition occurs when Q* 



1 and A — > 0. It corresponds to a growth rate 
a = ^/3fc(fc-l) 



(11.26) 



fe>0 



which vanishes. If a < 0, Q* = 1 whereas for a > 0, one has Q* < 1. 

The generating functions Gi and G2 can be introduced as in pi.3III.5P and the discussion of their long time 
behaviours is the same as in pi. 141) . For a flow u{t) such that dtu{t) = F{u{t)) and u(0) = uq, one can define an 
amphtude A as in (III.12P by : 



A(uq) = lim 



(11.27) 



The generating functions of the size Nt in the quasi-stationary regime are still given by pi.l6[) but the definition of 
A for continuous time leads to the following continuous time version of pi.l7p : 



As in (|lLT8l) . the expansions of A and F near the fixed point Q* can be related through 



A ' 

3 



2 V A 



F(3) 

~2A"' 



(11.28) 

(II.29a) 
(11.29b) 



The average size is given by (N)qs = —F^^"^ /A and diverges as one approaches the transition a — s- 0. For models with 
finite moments A;", the first moments (A^")qs obtained from the expansion pi.29p of A go to the same values as in 
(IIL22]) 

As pi.28[) is a first order differential equation, one can fully determine A using A'{Q'*) = 1 (which follows from 
definition (|lL28)) ') : 



AiQ) = (Q-Q:)exp 



A 



1 



Fiy) y - Ql 



dy 



(11.30) 



C. Universality of the stationary regime 

Near the transition, Q* is close to 1, A is small and for Q close to Ql, one can approximate F{Q) by the first two 
terms of its quadratic expansion near Q* (when /c^ <; go) : 

F{Q)^XiQ-Q:)+B{Q-Q:f. (11.31) 

Then pi.30p becomes near Q* : 



and from pi.l6|) . one gets that for /i small 



{^-'"Us — 2 . (11.33) 

(1 - 

Close to the transition, B ~ F"(l)/2 ~ (fc^ — k)/2 and this shows that in the quasi-stationary regime : 

{NU 0^ ^ (11.34) 

and the distribution of the ratio x — N/ {N)qs becomes universal [2^ : 

P{x) ~ 4x6-2^ (11.35) 

with the ratios of the first moments given by pi.22|) . 

The case of slowly decreasing branching rates (3k for which the second moment k"^ is infinite leads to other universal 
distributions near the transition, as shown in 40] and in appendix TK] 
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III. QUASI-STATIONARY REGIME FOR A BRANCHING RANDOM WALK IN PRESENCE OF AN 

ABSORBING WALL 

A. Evolution of the extinction probability 

We now consider a population for which each individual is characterized by its position x relative to an absorbing 
wall. Each individual diffuses and branches into k offspring with probability (3kdt during dt (for simplicity we set 
(3o = 0, i.e. there is no spontaneous death ; also there is no dependence on /3i as it has no effect on the population). 
The absorbing wall moves at constant velocity v and every individual crossing the wall disappears instantaneously 
(see figure [T]). In absence of the wall, the population grows exponentially and occupies a region which spreads in space 
linearly at a constant velocity Vc 0) HH . The presence of the wall introduces a competition between the exp onential 
growth and the absorption by the wall. For v > Vc, the wall moves faster than the population spreads |12l | and the 
population dies with probability 1. 

If one starts at time t — with a single individual at distance x from the moving wall, the extinction probability 
Qe{x, t) {i.e. the probability that all the descendants have been absorbed by the wall during the interval [0, t]) evolves 
according to : 

Qe{x + vdt^t + dt) = / — j^^Qe{x + rjv dt,t)dr] 



(III.1) 

J20k{Qe{x,t)'' -Qeix,t))dt. 



k=2 



To derive this equation, we divide, as in section |TT1 the time interval [0,t + dt] into a first time step dt and a second 
interval [dt^t + dt]. During the first infinitesimal time step dt, the individual diffuses (its position is shifted by an 
amount r]^/dt where 77 is a Gaussian random variable such that (77) — and {rf) = 2) and may branch into k > 2 
individuals with rate (3k- If this branching event happens, then, during the second time interval [dt,t + dt], all the k 
offspring present at time dt have independent histories. 



In the limit — > 0, this becomes 



dtQe = HQe) (ni.2) 



where the functional J- is defined as : 

J^iQ) = dlQix) - vd^Qix) + giQix)) (III.3) 

00 

fc=2 

Qe evolves according to a travelling wave equation [i^ . Here however the evolution is limited to a semi-infinite line 
with the boundary condition Qe{0,t) = 1 on the wall. The non-linear function g{q) contains all the information on 
the branching rates f3k of the walk. In the long time limit t 00, the extinction probability Qe{x, t) converges to the 
stable solution Q* of ^{Ql ) — and therefore satisfies : 

d^,Q:-vd^Q:+g{Q:{x))^Q. (III.5) 

For V > Vc, the wall moves faster than the spreading velocity Vc of the population, so that all the population gets 
absorbed and for large t, Qe{x,t) has the shape of a front moving at velocity v — Vc. Therefore Qe{x,t) converges to 
the unique stable fixed point Q* = 1 (the population becomes extinct with probability 1 whatever the initial distance 
X to the wall is). 

For V < Vc, there exists a non trivial fixed point Q* with Ql{0) — 1 and Ql{x) for a; — > 00 and Qe{x, t) — + Q*c{x) 
as i — !■ cx), with relaxation times : 

Qc{x, t) ~ Ql + Ai0i(a;)e-*/"i -|- . . . (III.6) 

These relaxation times t„ are related through A„ — — l/r„ to the eigenvalues A„ of the linear operator C obtained 
by linearizing the functional T{Q*c + e(j)) ~ e>C[0] + o(e) around the fixed point Q*. The eigenvectors (j)n{x) and the 
eigenvalues A„ satisfy : 



£[0„] = A„0„ 



(III.7) 
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with the linear operator C given by : 



with a potential V given by : 



- dliPix) + V{x)ii}{x) = (-A)V(a;) 
V{x) = '^-Y,h{kQl{xf-'-l) 



(IIL8) 
(IIL9) 

(III.IO) 



Defining ^'(s;) = (j){x)e ^^'^ maps ()IIL8[) to a Schrodinger equation 



k 

(note that for v > Vc^ one has Ql{x) = 1, so the potential V is constant and the spectrum is continuous. For w > Wc, 
Qe{x,t) converges to Ql{x) in a non uniform way, so that there are always some x for which Qe{x,t) — Ql{x) is not 
small and the long time behaviour of Qe{x,t) cannot be understood by hnearizing around the fixed point Q*f.{x) = 1). 

As V approaches Vc, the region where (5*(a;) is close to 1 increases (and diverges when v = Vc)- As a consequence, as 
V ^ Vc, the potential well V{x) in pil.lOp has a finite depth but an increasing length and this gives rise to a discrete 
spectrum for small eigenvalues. We will calculate the eigenvalues A„ of C in this limit v ^ Vc ui section IV Al and 
appendix |B] 

In the rest of this paper we will only consider the case v < Vc when the smallest eigenvalues of C form a discrete 
spectrum (only the construction of section HV Al will be also valid for v > Vc as one does not need there the slowest 
eigenvalue to be isolated). 

In constrast to (jll.3[) . the population is now not only characterized by its size Nt but also by the positions of all 
the individuals. One can however follow the same approach as in section |TT1 Let us introduce the generating function 
Gi : 

Nt 

Gi(x,t;/) = (ne-^^"-"^) (ni.ll) 

where / is now a positive test function of x and plays the role of fi in pi.3p , Nt is the size of the population at time t 
generated by an initial individual at x and the xf "' are the positions relatively to the wall of all its offspring at t. The 
independence of the lineages implies that Gi for the initial individual is the product of the Gi 's of its offspring after 
the first time step dt. Therefore, the evolution of Gi(x,t;/) can be derived as in (jlll.ip for Qe{x,t) and one gets : 

dtGi{x,t;f)=T{Giix,t;f)). (III.12) 

Absorption by the wall implies Gi{0,t; f) — 1 for alH > and, for a single initial individual at x when t = 0, the 
initial condition for Gi is : 

Gi(a:,0;/) = e--^(^). (III.13) 



B. Conditioning on survival 

As in section ini one can then try to calculate the size Nt of the population or the density p{x) at time t conditioned 
on the survival of m individuals at later time T = t + 1' (see fig. [2|) . To do so we introduce the two-time generating 
f miction at G2{x,t,t'; fi,^) defined as in (|II.5[) : 

Nt 

G2(x, t, i'; /, i/) = ( exp ( - ^ f{xf^) - lyNt+f) ) (III.14) 

i=l 

given that at i = there is a single individual at distance x from the wall. Once more, one can show that G2 {x, t, t'; /, 
(as a function of x and t only) evolves as Gi{x,t; /), i.e. according to piI.12p with the following initial condition : 



G2ix,0,t';f,iy)^e~f^^^Gi{x,t';i^). 



(III.15) 
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As in pi.8III.9III.10p . one can expand G2 in powers of e ^ : 

oc 

G2{x,t,t'-S,v) = Y,e--"'R^{x,t,t'-f) (III.16) 

m— 

00 

G2(x,i,t';0,i/) = Gi{x,t + t'\v) ^^e-'''^Pm{x,t + t') (III.17) 

m=0 

where Pm{x, t + 1') is the probabihty that the initial individual located at x has exactly m living descendants at t + 
and the generating function of / conditioned on observing a size m at t + t' is given by : 

<^exp ^/(a;f')) size m at t + i'^ = t, i'; /) (III. 18) 

where Rm{x, t, t'; f) is defined as : 

This shows that all the information about the regime < t < T, conditioned on having the final size Nt — m, is 
contained in the functions Rjn{x,t,t'; /). Eq. piI.15IIII.16IIII.17p are the analogues of pi.7III.8III.10p in sectionlTTl 

Inserting the expansions piI.16IIII.17p into piI.12p leads to differential equations for the i?m's and the Pm's. At 
lowest order, Rq and Pq satisfy piI.12p and, at first order, i?i and Pi are related to Rq and Pq through : 

00 

dtRi = dlRi - vd^Ri + ^ hikR^^^ - l)Ri (III.20a) 

k=2 

oc 

dtPi = dlPi - vd^Pi + J2 MkPt^ - l)Pi (III.20b) 

fc=2 

i?i(a;,0,t';/) = e'-^(^)Pi(x,t'), Pi(a;,0) = l (III.20c) 

By expanding further G2 in powers of e~'^ , one could obtain in the same way the evolution equations of the Rm's and 
Pm's and from pil.isp their initial values. 

C. The quasi-stationary regime 

In this subsection, we show how one can generalize the expression pL16p of section [IT] relating the properties of 
the quasi-stationary regime to those of the functional defined in pil.3p near the fixed point Q* for v < Vc- An 
alternative (easier) approach to calculate these properties will be described and used in sections |IV] and El 

The generating functions Gi and G2 evolve according to piL12p and, as soon as the test function / is positive, Gi 
and G2 converge to the non-trivial fixed point Q* of J-^ when v < Vc- In order to analyze the long time behaviour of 
Gi and G2, one can expand J-' around Q* : 

J'iQl + ec^) - T{Ql) + eC[4>] + ^.F'^) [0, -|- . . . (III.21) 



where the n-linear symmetrical functionals T^"-^ are defined for any arbitrary set of n functions V'i as : 

(III.22) 



T"- ^ Vl,---,V'n = 1 ^ 

dsi . . . dSn 



As discussed after (IIII.lOp . the linear operator L does not always have a discrete spectrum (in particular for v > Vc). 
We will assume here that v < Vc and that at least the smallest eigenvalue Ai is isolated. If the first eigenvalue is not 
isolated, all the construction below breaks down and the quasi-stationary state might not exist [l^ . 

As Qe, Gi and G2 all evolve according to piI.12p . one can study the dynamical system : 

idtu{x,t) ^T{u{x,t)) (III 23) 

\u{x,0) =uo{x) 
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When u{x,t) Ql{x) as t ^ oo and if C has a discrete spectrum for smaU eigenvalues, the convergence of u{x,t) is 
exponentiaL Then one can introduce as in pi.lSp the amphtude A of the first eigenvector 0i of C (with the largest 
relaxation time ri) : 

u{x, t)^Q*+ 0i(x)e~*/^M (uo(-)) + . . . (III.24) 

but now the argument of the amplitude A is a function (the initial condition ^0(2;) = u{x,0)). 

From the initial conditions piI.13IIII.15p of Gi and G2, the generating functions Gi{x,t;f) and 6*2(2;, t, i'; /, i/) 
thus have the following behaviours as t ^ 00 : 

Gi{x,t;f) = Q:(x)+0i(x)e-*/^M(e--^')+o(e"*/"O; 
G2{x,t,t'-f,y) - Ql{x)+<Pi{x)e-"^'A{e-fGi{-,t'-v))+o{e~"^'): 
Thus for large t and t', at leading orders, G2{x, t, t'\ /, v) becomes : 

G2ix,t,t';f,,^) ~ Q:(a;) + 0i(a;)e-*/^M(e-/(Q: + e-*'/^i0i^(e-'^))) 

~ Q:{x) + Mx)e-'^^'A{e-fQl) 



+0i(x)e-(*+*')/"^ —A (e-fQl + se-f^i) Aie-"") + ... 
as s=o 

As G2{x, t, t'; 0, 1/) = Gi{x, t + t'\ v), one can proceed as in pi.isp by writing A{e^'^) — AmC" 
and defined in piI.16IIII.17p for to > 1 : 



R^{x,t,t'-f) ~ A„(/.i(x)e-(*+*')/^i-^(e-/Q: + se-/0i) 

P™(x,t + t') ~ A„(/.i(x)e-(*+*')/^i 
Thus, the generating function of / given by piI.16p . when i' — > 00, has a finite limit 



s=0 



lim 

t,t'— ►CO 



Nt 



given by : 



exp(-^/(xf) 



= -A{e-fQ: 
qs as 







s=0 



(III.25) 

and get for Rm 

(III.26a) 
(III.26b) 

(III.27) 
(III.28) 



(note that as in pil.24p the argument of ^ is a function). We see that, as in pi.l6p . the properties of the quasi- 
stationary state are determined in pII.28P by the amplitude A (expression pII.28P does not depend anymore on t, t' , 
nor on the value of m used to condition the population size at t + 1' , nor on the position x of the initial individual). 
From the definition piI.23IIII.24p , one can see that the amplitude A satisfies : 



XiA{u) 



d 



Aiu + rTiu)) 



T=0 



(III.29) 



As in section ini this equation can be used to relate A and T. Note that equations similar to pil.29p appear in other 
contexts, in particular in the renormalization group theory where A is called a non-linear scaling field ^44i] . As in 
(IIII.2ip . one can expand A around Q* : 

A{Q: + ecj)) = AiQl) + eA^^^ [<^] + ^^'^^ [0, 0] + • ■ ■ (III.30) 

Successive derivatives of pIL29P as defined in pll.22p in the directions of the eigenvectors 0„ of C allow one to relate 
the expansions piI.30IIII.2ip of A and : 



A^'^[cj),,cj,,,<j>k] = 



1 



Ai — Xi — Xj 
1 

Ai — Ai — Xj — Afe 



(III.31a) 
(III.31b) 

(III.31c) 
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Then, from pil.28p . one can obtain the properties of the quasi-stationary regime. For example, the moments of the 
number of individuals in the quasi-stationary state are obtained by taking / = /i constant and taking successive 
derivatives at fj, — : 

asdfi M=s=o 
= A^^^ [0i] -I- ^(2) [^^ ,Ql] = l+ ^(2) , Ql] (III.32) 

= yt(i) [01] + 2y^(2) [01^ Q*] + _4(3) [0^^ Q*^ Q*] (ni.33) 

These expressions are exact and valid as long as the functional J-' has a non-trivial fixed point Q* and an isolated 
largest relaxation time ti. In order to go further, one needs to know more precisely Q* and 01. This will be done in 
section fVl in the scaling regime v Vc- 

IV. A MODIFIED PROCESS TO DESCRIBE CONDITIONED HISTORIES 

Before analyzing this scaling regime we are going to show that one can construct a modified process reproducing the 
history < f < T of the branching random walk of section Hill conditioned on a final size = 1- This construction 
is valid both below and above the critical velocity Vc- The quasi-stationary state for v < Vc discussed in section ITlI CI 
will appear as a particular case and its average profile will be calculated in section IIVBI 

A. Construction of the modified process 

One way of describing branching random walks conditioned on the survival of a single individual is to distinguish 
the path of the surviving particle from the other ones. In the mathematical literature [HI, US], this particle is called 
the spine and has its dynamics modified in order to prevent it from dying. The rest of the population is generated 
from this special particle. We show in the present section that this construction is general for branching random 
walks conditioned on having exactly one survivor at finite time T and requires only the knowledge of the extinction 
probability Qe{x,t) and the probability Pi(x,t) that there is exactly one survivor at t. 

To condition the evolution of the population on the survival of one individual at time t + t', equations piI.18l[TlI.19[) 
show that one should consider the ratio Ri{x,t,t'; f) = Ri{x,t,t'; f)/ Pi{x,t + t'). From (|III.12IIII.20aP satisfied by 
Rq and i?i, one can show that Rq and Ri satisfy : 



k=2 



d^R, = dlR, + [-v + 2^^^%ii-^~) d^Rr (IV.2) 



fe=2 



Pi{x,t + t' 

khQe{x,t + t'f-^ (r^-'Ri-r, 



(where we have used the fact that Qe{x,t) ~ Po{x,t)). The initial conditions (|III.20cP give Ro{x,0,t' ; f) ~ 
i?i(x,0,i';/) = e-^(-). 

We are now going to show that the functions Rq and Ri can be interpreted as the generating functions at time t 
of a modified process defined on [0, T] with T — t + 1'. To do so, we consider the following process for < i < T in 
the frame of the wall : 

• the system consists of a single particle of type Ai (the "spine" particle) and an arbitrary number of particles Aq; 

• a particle of type Aq diffuses at time t in the frame of the wall with a drift 
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and branches into k particles of type Aq at rate I3^\x, t, T) = PkQei^, T — t)'^ ^ : 

A, i!^^ kA, (IV.4) 



• the particle of type Ai diffuses in the frame of the wall with a drift 

v,{x,t,T)^-v + 2 p^^^ j,_^^ (IV.5) 
and branches into one particle of type Ai and k— 1 particles of type Aq at rate 13^'^ {x, t, T) — kl3kQe{x, T—t)''~^ : 

A, i!lJf:i^ A, + ik~ l)Ao (IV.6) 



One can notice that now the drifts and the mutation rates depend both on space and time through the functions 
Qe{x,t) and Pi{x,t). Moreover, as Pi(0,i) = and dxPi{0,t) ^ 0, the particle Ai is never absorbed by the wall 
(whenever it approaches the wall, it is pushed away (jlV.Sp from it by the drift dxPi/Pi which diverges when x ^ 0). 
The extinction probability Q'f\x,t;t') of a particle Aq is given by Q^^\x,t;t') ~ Qe{x,t) / Qe{x,t + t') and one can 
verify that it satisfies pV.l|l : indeed it is equal to Rq{x, t, t'; oo). Thus, at t' = 0, i.e. at the final time T, all particles 
^0 have disappeared. 

For this modified process, one can consider the new generating function G'"-' (x, t, T, /) = (exp(— /(s;^*'))) (resp. 
G^^-* (x, t, T, /)) where the summation at time t is over both particles Aq and Ai, given that we start with a single initial 
individual of type at position x (resp. Ai) at t = 0. In this modified process, the offspring are still independent 
as in (jIII.l[) and the function (resp. &^^) is the product of the G^'^'s of the offspring of the initial individual of 
type ^0 (resp. Ai). As in (jllLip . one can spht the time interval [0, T + dt] into two intervals [0, dt] and [dt, T + dt] 
and one gets : 

&'^'>{x,t + dt,T + dt,f) = 



^G(°) {x + r]Vdi + vq{x, 0, T + dt)dt, t, T, f)dr] 
PkQeix, T + dtf-^ (g(°) (x, i, T, ff - G(°) (x, t, T, f)] 



and, as dt is infinitesimal, G^"-* has the same evolution pV.ip as Rq. Similarly, writing the evolution of G'^-* in the 
same way shows that it has the same evolution as Since the G^^-* and Rq (resp. G^^^ and Ri) have also the same 
initial conditions, one has RQ{x,t,t'-J) = G^°\x,t,T,f) and Ri{x,t,t' ; f) = G'-^'^ {x,t,T, f) and the i?i's introduced 
in (jIII.19[) can be interpreted as the generating functions of a modified process defined on [0,T]. 

When T — > oo and v < Vc, the extinction probability Qe{x,T — t) converges to Ql{x) ^ and Pi{x,T — t) decreases 
as in piL26p . Thus, as T — s- oo, the drifts and the branching rates of the modified process become independent of t : 



v,{x,t,T) ^ wq{x) = -v + 2^^^^, f3l°\x,t,T) PkQlix)'-' (IV.8a) 



v,{x,t,T) ^ wiix)^-v + 2^^^^, p^^\x,t,T) ^ kPuQlixf-^ (IV.8b) 

(j)i{x) 

From the expressions pV.Sp . one can understand the best strategy for the system to have a single survivor at 
T = t + t' . The branching rates decrease as x increases and vanish as Ql{x) goes to : no population can develop 
in the region where Ql{x) ~ and it prevents the population from growing exponentially (which would not be 
compatible with a finite size at T). The particle Ai cannot be absorbed because vi{x,t,T) — s- oo as a; — s- 0. 

Remark : for the birth-death process of section III Bl one could similarly construct a modified process [2l| with 
branching rates PkQl^~^ for particles Aq and kPkQt'^^^ for the spine Ai, to describe the quasi-stationary regime. 

Remark : The construction of a modified process can be easily adapted to intermediate regimes < t < T 
conditioned on the survival of two (or more generally to p) survivors at T. The differential equations satisfied by R2 
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and P2 coud in the same way be interpreted as the dynamics of a system of three types of particles with the following 
branching rates : 

Aq kAo, rate /3fc(5e(x,r - t)''^^ 
Ai^ Ai + {k- l)Ao, rate kl3kQe{x, T - t)^-^ 



^ A2 + (fc- l)Ao, YaXekl3kQe{x,T -t) 



fc-i 



A2 ^ Ai + Ai, rate ^^^^ PkQe{x, T - t f-''Pi{x, T - t) 

where Aq and Ai have the same drifts as before and A2 has an additional drift 2dxP2{x,T — t)/P2{x,T — ~ 
'^dx'piix) / (j)i{x) for r ^ 00. The exponential decay of Pi to with a relaxation time ri shows that the reaction 
A2 ^ Ai + Ai occurs only for t ~T. 



B. Average profiles 

The average density profile p{X, t) of the population is defined such that p{X, t)dX is the average number of 
individuals located in the interval [X, X + dX] in the moving frame of the wall. It is easy to see that it satisfies : 

dtp^dlp + vdxP + Y.kh{.k-l)p (IV.9) 

with p(0,<) — because of the absorption by the wall and p{X,Q) = 5{X — x) where x is the position of the initial 
individual. If one introduces the growth rate a = f3k{k — 1), the solution is given by : 



{x-xy/4t _ -(x+xy/4t 



Similar expressions were obtained in [1] for a model of cell proliferation. 

If w > Vc = 2a^/^, the average density decreases to zero and this corresponds to an almost sure extinction of the 
population. For v < Vc, the divergence corresponds to the exponential growth whenever the population survives. 

We are now going to calculate the average profile PqsiX) in the quasi-stationary regime using the modified process 
of section IIV AI There are two contributions to this profile : 

Pqs(^) = Pl,st{X) + PoAx) (iv.u) 

where pi^st(-''^) and po^st{X) are the stationary average density profiles of the particles Ai and Aq in the modified 
process. In the moving frame of the wall and in the stationary regime, particles of type Ai (resp. Aq) have a drift 
wi{X) — —V + 2dx4>i/4>i (resp. wq{X) — ~v + 2dxQl/Ql)- The average density profiles po and pi satisfy equations 
similar to pV.9[) where we use now the drifts and the branching rates pV.8l) of the modified process : 

dtPo{X,t) = d],poiX,t)-dxiwo{X)po{X,t)) 

+Y,{k~i)hQ:{xf-'p^{x,t) 

k 

+ Y,k{k- l)l3uQl{Xf-^pi{X, t) (IV.12) 

k 

dtPi{X,t) = dj,pi{X,t) - dx {wi{X)pi{X,t)) (IV.13) 

There is no source term for pi since the number of particles Ai is conserved. As there is initially a single particle Ai, 
the density pi is the distribution of its position at time t. The stationary profile pi^st is directly obtained from the 
expression of wi{X) and PV.13P : 

pi,st(X) -C0i(X)2e-"^ (IV.U) 

where C is a normalization constant such that pi^Bt{X)dX = 1. The density pi.st depends only on the slowest 
eigenvector (pi and its shape near the critical velocity will be studied in section fVAl 

The extinction probability of Aq particles in the modified process is given by construction by Ro{x,t,t';oo) = 
Qe(x, t) / Qe(x, t + 1') and thus becomes Qe{x, t) / Q* {x) as ^' — > (X) in the quasi-stationary regime. Therefore a particle 
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Aq gets extinct with probability 1 in the long time limit. Thus particles Aq are produced by the particle Ai and later 
on are absorbed by the wall. Their stationary density po,qs satisfies : 

dj,po,st{X)-dx K(^)po,st(^)) /3feQ:(^)'~'(fc- l)Po,st(^) (IV.15) 

k 

= -g"{Q:{X))Q:{X)p^,,,{X). 
The change of variable po^st{X) — e^^-^ Q%{X)^j;{X) leads for V to an equation of the type : 

C[il)] = K{X) (IV.16) 
K{X) - -.g"(g:(X))pi,,t(X)e'^^ (IV.17) 

where C is the linear operator defined in (jlll.Sp . This equation is an inhomogeneous second order linear differential 
equation that can be solved easily (because one solution of the homogeneous equation, dxQl{x), is known). The 
general solution of (jIV.16p is 



dy ry 



[d.Ql{vWe--y A 



d^Ql{z)e-'"K{z)dz. 



The stationary profile po,st(-''^) is then obtained by choosing a = and 6 ^ oo for pQ^stiX) to vanish at X = and 
X ^ oo : ' 

PoAX) =e-"^Q:(X)a,Q:(X) 

dv r°° „ (IV. 18) 



V 



Equations pV.lli IIV.141 IIV.18|) give the exact average quasi-stationary profile in terms of only Q* and which 
have been used for the construction of the quasi-stationary state. Section |V] will describe the limit v Vc of Q* (x) , 
(x) and pqs j X). 

The results (lIV.lll IIV.141 IIV.18P could also be derived from the dynamical system approach of section IIII CI The 
average quasi-stationary profile pqs{X) at distance X from the wall is obtained by considering the test function 
f{x) — pS{x — X) and taking the /i-derivative of pil.27p a.t p = 0, as for the average size in pil.32p . Relations 
(IIII.31al IIII.31bp between A^^''[(j)i, 4>j], C and JF^^' (/j^J are then exactly equivalent to the differential 

equations (jIV.12|IV.13P in their stationary regime. 

V. UNIVERSALITY NEAR IN THE CRITICAL VELOCITY 

In this section, we determine the scaling properties of the quasi-stationary profile when Vc ~ v ^ 0^ ■ We are going 
to see that, in this limit, everything can be expressed in terms of the shape Qv^ix) of a front on the infinite line at 
the critical velocity. 

A. Extinction probability and relaxation times 

The stable fixed point Q* (x) of the functional defined in piI.Sp satisfies the differential equation pil.Sp with the 
boundary condition Ql{0) — 1 and Ql{x) goes to or 1 as a; ^ oo. 

In appendix [Bl we use a perturbation theory to analyze pil.5l) on the infinite line for Vc — v 0^ to obtain the 
extinction probability Q*{x) in presence of the wall in the same limit. The solution Qy^ of pil.sp on the infinite line 
such that Qv^{xo) — 1/2 (where xq is some arbitrary position), Qy^{x) ^ 1 as x — *■ —oo and QvS^) ^ as a; — *■ -|-c» 
is given for large negative x : 

Qy^{x) - 1 + {A,x + B,)e-^-/^ + O (e"=-) (V.l) 

where Ac and Be are constant depending on xq- For v close Vc, one can show (see appendix [B]) that, on a domain of 
length L that we will call region I, Ql{x) is given by : 



ArL . /irx 



region I 



(x) ~ 1 ^sin — exp -^{x - L - Bc/Ac) (V.2) 
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where Ac and Be are the constants defined as in (jV.ip . The length L is given by : 

L = - - {vc-v)-^/^ (V.3) 

V — >Dc 

and diverges at the critical velocity. In the domain x > L (region II), the non-linearities of (jlll.sp must be taken into 
account and Q* is given at leading order by : 

Qlix) = QvA^-L-Bc/Ac) + 0{l/L^) (V.4) 

region II 

The relaxation times r„ of the dynamics (jIII.2[) near the stable fixed point Q* are related to the eigenvalues 
A„ = — 1/ r„ of the linear operator piI.Sp . An eigenvector (pn with eigenvalue A„ satisfies the second order differential 
equation : 

9|(/)„ - Vdx(f>n + g' {Qvix)) cj)n = K4>n (V.5) 

with the boundary conditions : 

0n(O) = 0, 4>^{x) (V.6) 

As for Q* , the shape of the eigenvectors (pn can be obtained from the ones on the infinite line as shown in appendix 
(|B|) . The effect of the boundary condition (|V.6p on the wall is to select the ones which vanish at a; = 0. A perturbative 
expansion in Vc ~ v and A shows that the first eigenvalues A„ in presence of the absorbing wall are given by : 

One can notice that the first three terms in this expansion are independent of the precise form g{Q) of the non- 
linearities. The eigenvector 0„ associated to A„ is given (up to a multiplicative constant) at leading order by : 



4>n[x) ~ sm(— — exp — [x- L- — \ 

region I 2n7r \ L J L^V A^ J 



(V.8a) 



in the region I of size L near the wall and by : 



0„(x) ~ d.,Q,Sx- L- BjAc) (V.8b) 

region II 

in the region II such that x> L. Appendix IB] also gives the order of magnitude of the first correction to (pn in both 
regimes x <^ L and x > L. As the relaxation times are given by r„ = — 1/A„, the longest relaxation time ri is given 
dVJl by 

n^^^^(.,-«)-3/2 (V.9) 



47r2 y/2v^ 

and is much larger than all the other r„'s {n > 2). One can notice that the same expression of ti was obtained in [l^ 
by a very different approach. 

B. Quasi-stationary profile 

The quasi-stationary profile pV.lip is the sum of two contributions (IIV.14IIV.18p of the modified process. For v 
close to Vc, the normalization constant C in ()IV.14p is dominated by the contribution of region I and is given at lowest 
order by C ~ 87r^e"'^'^-^+^'=/^^V(^c^c-^^)- Thus the stationary profile pi^st(-''^) at distance X from the wall in the 
modified process is given by : 

pi^X) -sin — (V.IO) 

region I L, \ L J 



for X < L. For A" > L (region II), the density is given by : 



^ TT ^" \2 2r3 e"""" [^-Q-c {X-L- Bc/Acf (V.U) 

region II Aiv^L'^ 
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The corrections to this profile are of order \/L in both regions. 

The density profile /9o,st(^) given by pV.18[) can be rewritten in the following way 



X 



dy 



[d.Qt{yWe-v 



{loo - l{y)) 



where 



and la 



d,Ql{z)g"{Ql{z))pi,st{z)dz 



(V.12) 



(V.13) 



lim 



y^+oc I{y)- For X < L, 01 (x) and dxQe are of order e"^'^ ^-'/^ and Q is of order 1 ; for x — L ^ 1 



(region II), all the terms in the integral (jV.13[) decrease exponentially. Thus, the integral is dominated by the region 
of size of order 1 near a; ~ L and /no scales for Wf, — w — > O"*" as : 



1 

l3 



Sir' 



[d:,QM?9"{QvM)e-'"^'dz 



(V.14) 



One notices that I^o is negative since Q* (a;) is a decreasing function and that loo is independent of the choice of the 
reference xq such that Qv^ixo) = 1/2. On the other hand, for the same reasons, I{y) is exponentially small and can 
be neglected compared to I^o in the linear region < y < L. 
At leading order in the region I, one gets from (|V.12p : 



X 



[d.Q*M?e-y 



dy 



(V.15) 



From (|V.2|) . dxQl{x) is given for a; < L by 



dxQtix) 



ArL [V. 



. /'!TX\ 

— sm — 
2 \ L J 



TT /■KX\ 

— cos 

L \ L ) 



region I TT 

up to exponentially small corrections. The integral in (jV.15[) is equal to 



X 



1 



[dxQl{y)fe-v 



dy 



^2^v,(L+BJA,) j-X 
^v^(L+BJA^) 



dy 



[i^sin(^) + fcos(Y)]' 
sin (2^) 



Al fsin(^) 



Finally the density of Aq particles is given in the region I by : 



yOo,st(^) 



ArTT 
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exp 



rcos(^: 



(V.16) 



In the region x < L, the number of particles Aq is exponentially large and pi^st(^) is bounded and much smaller than 
Pa.stiX) as shown in (jV.lOl |V.16[) . Thus the quasi-stationary average profile pqs{X) scales in the linear domain as : 



^e--(^-^)/^sin 



region 



I 4ttL^ 



TrX\ 



with a constant K given by 



K 



AM 



dxQvSz)Y9"{QvSz))e-''''dz 



(V.17) 



(V.18) 



One can check that this expression of K is independent of the reference point xq chosen for Q . For X ~ L oi order 
1, pqs(X) is of order 1/ Li^ and in the domain X > L (region II), the density decreases exponentially, so that the total 
number of individuals is dominated by the region < X < L. The average size in the quasi-stationary regime is thus 
given by : 



{N)c 



Pqs{X)dX 
{vc - vf^"^ exp 



K 



(V.19) 
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The divergence (IV.19P of {N)qs can also be obtained from (|III.32|) by expanding Q* in the basis of the eigenvectors 
and using (|III.31al IIII.31b[) : the leading term (|V.19p could then be obtained by truncating the expansion of Q* after 
the first eigenvectors 0i, as the first eigenvalue Ai ~ {vc — u)^'^/^ is much smaller than the next ones that scale as 

In Ji, iSj, Brunet et al. considered a population whose size N is kept constant by always selecting the N individuals 
with the largest Xi. From they predicted that the effect of the finite size TV is to shift the average velocity (v) of 
the population by an amount : 



21og^iV ■ ' - ' log'N 



I \ ^'^^ , ( 2N31oglogiV 



This TV dependence of the velocity was recently derived rigorously in [33| . It is interesting to notice that this relation 
between (w) — Vc and iV is the same as (jV.19[) between (A^) and v — Vc 'va. the quasi-stationary state of the present 
model with an absorbing wall. 



VI. EXPONENTIAL MODEL 



There is a simplified version of the problem of a branching random walk with an absorbing wall, the exponential 
model, which can be solved exactly [3 0]. In this exponential model, time is discrete and at each generation, each 
individual is replaced by its offspring distributed according to a Poisson point process : if the parent is at position x, 
then for every interval [y, y + dy\ there is an offspring in this interval with probability '\\}{y — x)dy. 

By analyzing the first time step one can show that the extinction probability Qe{x, t) at time t of the descendance 
of an individual located at distance x from the wall at time t — evolves according to : 



(VI.l) 



Qe{x, t + l)= T{Q^{x, t)) = exp '4>{y + v-x){l- Qe{y, t)) dy^ 

For a general -0, is a functional of Qe for which one could try to generalize the approaches of sections IlIII and [TVl 
For the special case = e"^, however, the problem reduces to the study of a one-dimensional mapping. For any 
Qe{x, 0), it is easy to see from (jVI.ip that a,t t > 1 Qe(x,t) takes the form 

Qeix,t) ^ e-'^^'' (VI.2) 

and depends on a single parameter Cj which satisfies the following recursion : 

ct+i = F{ct) (VI.3) 

F{c) = ( 1 - / e-ye-'^^'dy) (VIA) 



For all V, there is an attractive fixed point > and the population always has a non-zero survival probability 
1 — Ql{x) = 1 — e"'^*'^ in the long time limit. The relaxation of Qe{x,t) to Ql{x) is entirely controlled by the 
relaxation of Ct to c* : there is a unique eigenvalue A = F'(c,) < 1 with the eigenvector 

(j)i{x) ^ - exp{x - 0^,6"^) (VI. 5) 

and one has Qe{x, t) ~ Q*e{x) — KK^e^^'^''^ . Generating functions Gi and G2 defined as in (jIII.lllIII.14p also satisfy 
the same dynamics (IVI.ip . A general initial condition Q{x, 0) = Qo{x) is mapped at the first type step to a function 
Q{x^ 1) — e^^^'^o)'^ where the coefficient c((5o) is defined by 



c(go) = e-" (^1 - e-yQo{y)dy^ (VI.6) 

and has the following long time behaviour : 

Q{x, t) = exp(-cte^) ~ Ql{x) ~ S(c(go))A*-ie-=*^^ (VI.7) 
with an amplitude B{c) similar to the function A of section [III As in pi.l7p . the function B{c) satisfies : 

A6(c) = 6(F(c)), B'{c,) = 1 (VI.8) 
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with F given in (|VI.4p . 

The reduction to a one-parameter family exp(— ce^) makes the analysis of this spatial model similar to what we did 
in section nil 

For example, the generating function of the number N{X) of individuals at a distance x greater than a given X from 
the wall is given in the quasi-stationary regime by (jIII.28|) with A{Q) = B{c{Q))/ A and f{x) — fj,9x{x) = iiO{x — X). 
For the exponential model, it can be reduced to a simple function of /i : 

By decomposing c ((QJ + s<j)i)e~^^^^ in the following way : 

+e'\l ~ e-^) r e-y {QHy) + scp,{y)) dy 
Jx 

= + sA + (1 - e-^) (Jo(^) + sAJiiXj) 
Jo(X) = e-" / e-y^'^'dy (VI.9) 



Jx 

MX) = - / e-^'^'dy (VI.IO) 

where one can verify that Ji(0) = 1, one obtains the generating function : 

(e-'^^(^))q, - (1 - (1 - e-'^) Ji(X))6' (c, + (1 - e-^)MX)) (VI.U) 

It follows that the average density profile {p{X))^s obtained by taking the derivative of {N{X))^s and the average size 
of the population are given by : 

(p(X)>q, = -J[{X)+B^^\c^MX) (VI.12) 
(iV)q, - 1+ (-S(2)(c,)] Jo(0) (VI.13) 

Let us now analyze the large v limit (which for the exponential model plays the role ^ of the f — > Uc limit of 
section IVT) . The mapping F{c) has the following expansion near c = (corresponding to Q%{x) — 1) : 

F(c) = e-" (^-clnc + (1 - 7)c + ^ - ^ + . . }j (VI.14) 

where 7 = — r'(l) is Euler's constant. For large v, the stable fixed point c* of F is given by 

c* ~ e"^"^ exp(— e") 

and A ~ 1 — e"". As in section|lTl one can compute B^"^^ ~ — 1/c* from (|VI.8|) . The functions Jq and Ji in (|VI.9IVI.10p 
have different shapes for X < L (region I) and X > L (region II) where the length L is given by : 

L = -lnc, -e" (VI. 15) 

Thus the leading terms in (|VI.12IVI.13p when X < L (region I) correspond to Jq contributions and, for v — > 00, the 
average quasi-stationary profile and sizes are given by : 



1 



(p(X))q, ~ 1''- ' (region I, X<L) (VI.16) 

(A^)qs ^ —^je^ (VI.17) 

In region II (for X > L), {p{X))qs is decreasing as e""^*^ . As in section |Vl the region of length L near the wall has 
the dominant contributions to the size of the population. As in section IIVI one could interpret the Ji terms as the 
contribution of the Ai particle and the Jq terms as contributions of the Aq particles in the context of the modified 
process. 
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One can notice that, in the exponential model too, (jVI.17|) yields the same relation between the size of the population 
and the velocity as (39) in Q (where the size of the population is kept constant but the velocity fluctuates). 

In the limit v — > oo, it is possible to obtain the whole generating function (e~^^)qs. In this limit, N scales as 1/c* 
and one should take /i of order c* and thus B{c) needs to be known for c — c^, — 0(c»). If F{c) in (|VI.14[) is truncated 
after the first two leading terms near c, as in section fll CI and appendix [Xl one obtains yB(c, + c,it) ~ ln(l + u) in the 
limit V —^ oo and the following generating function of the size : 

/ qs ti^oo 1 + fj. 

Thus in the quasi-stationary regime the size of the population has an exponential distribution in contrast to what 
happens pi.33[) in the Galton- Watson process. In particular, the ratio {N'^)qs/ {N)"^^ goes to 2 : this result is similar 
to numerical results for lattice branching random walks obtained in p^ . 

VII. CONCLUSION 

In the present work, we have studied the quasistationary regime of a branching random walk in presence of an 
absorbing wall below the critical velocity [l2|. To do so, we have developed two methods. The first one, discussed 
in section IIIIl is a dynamical system approach allowing to determine the generating functions of the size of the 
population. 

In the second approach, one constructs from the original process a modified stochastic process equivalent to the 
quasi-stationary regime (section HVT) . This construction requires the knowledge of Po{x,t) = Qe{x,t) and Pi{x,t), 
the probabilities of observing and 1 survivor at time t for an initial individual at x. These methods allowed us to 
determine ()V.17|V.18IV.19|) the average population size and the density profile in the quasi-stationary regime near Vc- 
The average profile has a universal shape (|V.17|) in a domain of size L ^ {vc — w)~^/^ and the average size diverges as 
in (|V.19p in a universal way which does not depend on the precise form of the branching rates Pk- These quantities 
have also been obtained for the exponential model in section IVIl We noticed that both for the branching random 
walk and the exponential model, the relation between the velocity and the size of the population (jV.19IVI.17p is the 
same as in a model recently studied in [!|| for a population of constant size. 

Beyond the average quasistationary profile which has a universal shape near Vc (i.e. is independent of the precise 
form of the non-linearities g{q) in (jIII.4[) ). it would be interesting to see whether the whole quasi-stationary measure 
is universal when v — > Vc- In particular, it would be interesting to know whether the size of the population has an 
exponential distribution in the quasi-stationary regime as in the exponential model IjVI.lSp . 



Our approach was limited to the case v < Vc- For v > Vc one expects [12| that there is no quasi-stationary regime. 
The construction of the modified process of section IIVI remains however valid and this could be a starting point to 
understand the absence of a quasi-stationary regime. 

More generally, the construction of the modified process is valid as long as the individuals are independent. In 
particular, a modified process could be constructed to describe the evolution of a population conditioned to a finite 
size Nt = 1 in non-uniform media, as in [3ii i4i]i where the branching rates, the drift and the diffusion coefficients 
depend on the position x and where the domains have more complicated geometries (c?-dimensional with absorption 
on the boundaries). This could be a way of understanding the effect of the geometry of the domain on the existence 
and on the properties of a quasi-stationary regime. 

Using the modified process, one could also study the dynamical properties and the correlations in time of the 
quasi-stationary regime. In particular, one could try to determine the statistical properties of the genealogies in this 
quasi-stationary regime and compare them with the results obtained recently for a population of fixed size 0, ■ 

Appendix A: UNIVERSAL DISTRIBUTIONS OF POPULATION SIZES IN THE QUASISTATIONARY 
REGIME NEAR THE TRANSITION IN BIRTH-DEATH PROCESSES : THE CASE OF SLOWLY 

DECREASING BRANCHING RATES 

In this appendix, we discuss briefly the quasi-stationary regime of the Galton- Watson process when the branching 
rates Pk decay too slowly for k'^ to be finite. We consider the case where for large k 

(3k ^ Bk-^-^ (A.l) 

with 1 < ?7 < 2. Then, for Q close to 1, F{Q) takes the form for Q close to 1 : 

F(Q)~a(g-l) + C(l-Q)" + ... (A.2) 
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with C = Br(-?7) and a = - l)/5/c. 

For a < 0, the fixed point Q* = 1 is attractive and one gets from (|II.30p for Q close to 1 : 

A{Q) ~ ^^'^^ (A.3) 

(l_g(l_Q).-l)V(.-l) 

For small a > 0, the fixed point becomes 1 — Q* ~ (a/C)^/'-''^^-', the eigenvalue F'{Ql) = X scales as A ~ — — 1) 
and by integrating (|II.30|) . one gets : 

A{Q) ^ (1 - g)!-" (l - f (1 - Qr-') (g) • (A.4) 

Then by taking /i small in pi.l6p . one gets for a < : 



whereas for a > 0, one gets : 



(e-^^)q,^ + ) (A.6) 



Expressions (|A.5IA.6p . valid for \a\ small, show that both below and above the transition, the distributions of N in 
the quasi-stationary regime are universal (as up to a rescaling they depend only on rj) and that they differ from what 
was found (|II.33p in the case 77 > 2). 



Appendix B: PERTURBATIVE CALCULATION OF THE SHAPE OF THE FRONT, OF THE 
EIGENVALUES AND EIGENVECTORS FOR v CLOSE TO Vc 

In this appendix we determine perturbatively, for v — Vc small, the stable solution Q* of pil.5[) , the eigenvalues A„ 
and their associated eigenvectors of the linear operator £ defined in (jlll.Sp . 

For the shape of Q* as for the eigenvectors our main idea is to derive a perturbation theory in powers of Wc — w (and 
also A for the eigenvectors) in the region where the non-linear terms cannot be neglected (i.e. in the region where 
1 — Q* is not small) and to match this expansion with the solution which can be obtained in the linear region (where 
1-Q* < 1). 

The shape of Q* : 

It is convenient to consider the solution Qv{x) of pil.sp 

Q':,~vQ',+g{Q,)=0 (B.l) 

on the infinite line such that Qv(x) — > 1 as a; — )■ — c» and Qvix) — > as x — > +00 (instead of the solution Ql{x) on the 
semi infinite line). Here g{Q) = Xlfc PkiQ'' — Q)- The boundary conditions at ±(X) determine Qy up to a translation. 
One particular solution may be selected by imposing that 

Qv{xq) = i (B.2) 

where xq is a fixed position (which can be chosen arbitrarily and of course will play no role in our final results). For 
V < Vc — 2-\/(?'(l), the solution of (jB.l|l has damped oscillations, as a; — > —00, of the form 

Q,{x) = 1-1- f/„ sin f ^^^^±M\ exp -t- O {e^ (B.3) 



L J ' \ 2 
where the length L is defined by : 

L = 2tt{vI-v'^)-^ (B.4) 



21 



so that for large L, one gets 



27r2 



V ~Vc 



.i2 



(B.5) 



The constants Uy and $„ in (jB.3[) are a priori comphcated functions of v (or L) and also depend on a^o- 

If is the position of the right-most zero of 1 — (x) , then the solution Q* of (jIII.5|) we are looking for is given 
by: 



lix) = Qy{x + ul) 



(B.6) 



As V approaches Vc, the length L diverges, and (|B.3P allows one to estimate ul up to corrections exponentially small 
in L 



VL 



(B.7) 



We now assume that the solution Q^^ of (|B.1|1 is known for v ~ Vc- One can then expand Qv{x) in powers of or 
in powers of f c — w 



1 V 

Qv = Qv^ + .... = Qv^ + {vc-v)-^R + 



Putting (jB.Sp into (jB.l|) one gets that R should satisfy 



R" ~v,R' + g'{Q,^)R 



27r2 



(B.9) 



One can solve this equation with the required boundary conditions (the equations leading to higher order terms in 
the expansion could be solved as well) and one gets : 



R{x) = — g;^(x) 

Vr 



dy 



+ 00 



For a:; ^ — oo, one can show from (|B.1|) that 

Q,^{x) = 1 + {A,x + B,)e^ + H2(x)e''=" + H3{x)e^ + ... 



(B.IO) 



(B.ll) 



where the coefficients Ac and Be depend on xq and the polynomes could be obtained explicitly in terms of these 
coefficients Ac and Be by analyzing the non- linear terms of (jB.ip . 

One can also show from (jB.lOp or (even better) directly from (|B.9p that for x ^ — oo : 



R{x) = ttM - 



+ Cx + D\ e— +0(6"'=^) 



(B.12) 



where the constants Ac and Be are the same as in (|B.lip and the constants C and D could be determined from the 
knowledge of Q^^ [x) and (jB.lOp . 

This leads (|B.8IB.lllB.12p for large negative x to 



Acx + Bc + j-^ 



6 



—x^ + Cx + D 

Vc 



Qv{x) = 1 + 

+ 0(6^"=^) 

On the other hand for large negative x, in the range 1 ^ |x| ^ L, expression (|B.3p becomes 

"7r(x + $„) 7r3(x + $„)3 TT^x{x + ^y) 



}y{x) = 1 + Uy 



or 



L 6i3 VcL^ 

The comparison of (jB.lSp and (|B.14p leads to the following expansions for Uy and $^ 



+ 



(B.13) 



(B.14) 
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Then using (|B.6IB.7p one gets the following expressions for Q* (x) in the region where L — x 1 



Qlix) = 1 - exp 



and in the region x > L oi x — L = 0(1) 



v,{x~L~^) 



Trx\ 

t) 



o 



o 



X — L 



(B.16) 



(B.17) 



We emphasize that the knowledge of Q^^ {x) is sufficient to determine all the higher order corrections in expansions. 
The eigenvalues and the eigenvectors (j>\ : 

We consider now the shape of the eigenvectors (f)\,y of the linear operator C on the infinite line. On the infinite 
line, to each A, one can associate an eigenvector (l)x,v which satisfies : 

Kv - ^^'x,v + 9'{Qv)^x^v = (B.18) 
with the boundary conditions (|V.6|) 0a, 11(2;) — > as x — > zLoo. For v < Vc, the solution when x — > — cxd has the form 

7r(x + ^'A,i,) 



4'\,v{x) = Va,i, sin 
where the length La.d is defined by 



exp 



) +0(6"") 



Lx,y = 27r(< - 4A - v^y^ = L 1 - 



(B.19) 



(B.20) 



On the other hand, it is easy to see that Q'^ (x) is an eigenvector for A = and v — Vc, (in fact for any v, one has 
^0,1; (2;) — Qy{x) up to a multiplicative constant) and one can try to expand c/fx^y in powers of A and Vc — v 



A" 



L2" 



(B.21) 



n+m>l 

The Snm^s could be determined recursively solving inhomogenous equations similar to (|B.9p . For example 

Sia-VcS[o + g'iQvJSiQ = Q'^^ 
S'oi-VcS'o. + g'iQ.JSoi = Q'' ~ 9"{QvjRQ' 



As in (|B.13[) one gets that for large negative 



BrVr 



A, + OiX) + 0[ 



+ 0(e"=") 



(B.22) 



Comparing (|B.22p with (jB.19|l in the range of large negative x with x ^ Lx,v (as in (|B.13IB.14p ) one gets that 



Vx.. 



Be 
Ac 



O 



Vc 



AcVr „ 



O 

o 



r2 



Lx,v 



(B.23a) 
(B.23b) 



On the infinite line, as A varies, all the eigenvectors 4>x.v satisfy the boundary conditions at x ±00. In the 
presence of a wall, the eigenvector 0a, u has to vanish at the wall, therefore the rightmost zero x — —L — of Qv{x) 
must coincide with a zero x — —nLx^v — ^\,v of 4>x.v The boundary conditions select that way a discrete set of 
eigenvalues A„ . For Vc — v small, one gets then, up to exponentially small corrections in L : 



nL 



X,v 



(B.24) 
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where n is an integer larger or equal to 1. Using the expressions (jB.15IB.23a|) and (jB.20p this becomes 



n 1 - 



and this leads to 

One should notice that, for n = 1, Ai L"^ whereas all the other eigenvalues A„ scale as L~^. 

In the frame of the wall, the n-th eigenvector (/)„ is, up to exponentially small corrections, given by : 

(t)n{x) ~ 4'X„,v{x + Ul) 

and therefore in the range where L — x <C 1 one has 



(B.25) 



<l>n{x) 



2m: 



mTX\ „ / 1 



o 



and in the range where x — L is of order 1 



hn{x) ~Q'^^ [x- L~ ^ 
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